function x = sollu(A,b)

	L = [];
	U = []; 
	x = [];
	[L,U] = lu(A);
	y = soltrinf(L,b);
	x = soltrsup(U,y);
	x';
	 
endfunction
